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ABSTRACT 

We present the distributions for gluon radiation off stop-antistop particles produced in e + e~ 
annihilation: e + e _ — > tig. For high energies the splitting functions of the fragmentation 
processes t — > tg and g — > it are derived; they are universal and apply also to high-energy 
stop particles produced at hadron colliders. 



Introduction. Stop particles are exceptional among the supersymmetric partners of the 
standard-model fermions. Since the top quarks are heavy, the masses of the two stop particles 
t\ and t 2 ? mixtures of the left (L) and right (R) squarks, may split into two levels separated 
by a large gap |1|]-[[|. The mass of the lightest eigenstate t\ could be so low that the particle 
may eventually be accessible at the existing pp and even e + e~ storage rings. So far the result 
of search experiments at e + e~ colliders |4], [5| has been negative and a lower limit of 45.1 GeV 
has been set at LEP |5j for the L/R mixing angle outside the band of cos 2 9 t between 0.17 
and 0.44 and for a mass difference between the t± and the lightest neutralino x5 °f m ore 
than 5 GeV. The higher energy at LEP2 and dedicated efforts at the Tevatron will open the 
mass range beyond the current limits soon. 

To begin, we briefly summarize the well-known theoretical predictions for the cross section 
of the production process [Fig.|l](a)] 

e + e~ — > ti t\ 

For a given value 9t of the L/R mixing angle, the vertices of the t\ pair with the photon and 
the Z boson may be written as ie^Qlp^ ~ PtJ/^ where pi x and p% are the 4- momenta of the 
stop and antistop squarks, and the charges read 

Q-y = ~e t 

Qz = (cos 2 6 t — 2e t sin 2 8w) / sin 29w 

respectively. 9w is the standard electroweak mixing angle and eo = vivra is the electromag- 
netic coupling to be evaluated with a (Mz) = 129.1 in the improved Born approximation 
||. The Z boson coupling vanishes for the L/R mixing angle cos 2 9t — > 2etsm 2 9w ~ 0.30. 
Defining the 7 and Z vector /axial- vector charges of the electron, as usual, by e e = — 1, 
v e = — 1 + 4 sin 2 9 W and a e = —1, the cross section can be expressed in the compact form 

, 2 (v 2 e +a 2 e )Q z s 2 
7 4sin 2 2^ (s-M|) 2 + M|r| 

VeQ-yQz S{S-M 2 Z ) 

s\n29 w (s - M\) 2 + Mfr| 

where \fs is the center of mass energy and Mz-, r \ are the mass and the total width of the 
Z boson, respectively. The P-wave excitation near the threshold gives rise to the familiar 
P 3 suppression, where (3 — (1 — 4m? / s) 1 / 2 is the velocity of the stop particles. Angular 
momentum conservation enforces the sin 2 9 law, das/d cos 9 = | sin 2 9, for the angular 
distribution of the stop particles with respect to the beam axis. 

QCD corrections. Gluonic corrections modify the cross section ||q The virtual cor- 
rections, Fig.|T](b), can be expressed by the form factor 



1 Since we focus on QCD gluon effects for light stop particles in the LEP range, we do not take into 
account quark-gluino loop effects, assuming the gluino to be heavy; these loop effects have been discussed 
for squark production at the Tevatron in Ref.|9) and at e + e" colliders in Ref . []Tc| . 
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where a s is the strong coupling constant and the kinematical variable w is defined as w = 
(1 — f3) /{I + The form factor is infrared (IR) divergent. We have regularized this 
divergence by introducing a small parameter A for the gluon mass. The IR singularity is 
eliminated by adding the contribution of the soft gluon radiation [Fig.[I](c)], with the scaled 
gluon energy integrated up to a cut-off value e g = 2E g mt /^/s <C 1. The sum of the virtual 
correction (V) and the soft-gluon radiation (S) depends only on the physical energy cut-off 
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After including the hard gluon radiation, the dependence on the cut-off e g disappears 
from the total cross section. The total QCD corrections can finally be summarized in a 
universal factor 

e^e -> ti h (g)\ = a B 
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with (Fig.|) 
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Very close to the threshold the Coulombic gluon exchange between the slowly moving stop 
particles generates the universal Sommerfeld rescattering singularity [11] jT — > 7r 2 /2/5, which 
damps the threshold suppression, yet does not neutralize it entirely. Employing methods 
based on non-relativistic Green's functions, an adequate description of stop pair production 
near threshold has been given in Ref . |12[ , which also takes into account screening effects due 
to the finite decay width of the stop particles. In the high-energy limit || the correction 
factor in eq.(|3|) approaches the value (1 + 4a s /7i). 

In this note we present a general analysis of hard gluon radiation. We also include stop 
fragmentation due to collinear gluon emission in the perturbative regime at high energies 
and we give an account of non-perturbative fragmentation effects. 
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For unpolarized lepton beams the cross section for gluon radiation off t\ squarks 

e + e" -^iiixg 

depends on four variables: the polar angle 6 between the momentum of the t\ squark and 
the e~ momentum, the azimuthal angle \ between the tit±g plane and the plane spanned by 
the e ± beam axis with the t\ momentum [see Ref. and two of the scaled energies x(ti), 



x(t\), z(g) in units of the beam energy. The energies are related through x + x + z = 2 
and vary over the intervals fi < x,x < 1 and < z < 1 — /x 2 , where \x = 2m$ ,/ \fs denotes 
the squark mass in units of the beam energy. For the angles between the squark and gluon 
momenta we have 

2 — 2(x + x) + xx + a 2 
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The spin-1 helicity analysis of the cross section results in the following well-known angular 
decomposition [14| 

da 3,„ 9 ^ da u 3 . 9 . da L 

-(1 + cos 2 6)- — - + - sin 2 0- 



dxdxd cos 6 dx/2n 8 dxdx 4 dxdx 



da 1 3 . 2 <io" T , . 

sm2t'cosx3 — — + — sm t^cos2x- — — (4) 



2a/2 dx<ix 4 (ixdx 

[Z7 = transverse (no flip), L = longitudinal, / = trv*long, T = trv*trv (flip)]. If the polar 
and azimuthal angles are integrated out, the cross section is given by a = a u + a L . 

It is convenient to write the helicity cross sections as 

I3 3 da j a s S j + jj 2 N j 



as dxdx 4-7T (1 — x) (1 — x) 



(5) 



The densities and are summarized in Table |l|; p is the momentum of the t\ squark, p 
and k are the longitudinal momenta of t\ and g in the t\ direction, and px is the modulus of 
the transverse t\, g momentum with respect to this axis [all momenta in units of the beam 
energy]. Since /, T correspond to 7, Z helicity flips by 1 and 2 units, they are of order p T and 
Py, respectively. Note that the threshold suppression is absent in the U, I, T components 
and attenuated in the leading longitudinal L term as expected from eq.(f|). 



Fragmentation. In the limit where the gluons are emitted from fast moving squarks with 
small angles, the gluon radiation 

h -> tig 

can be interpreted as a perturbative fragmentation process. From the form of the differential 
cross section dajdzdp\ we find in this limit for the splitting functions, in analogy to the 
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Table 1: Coefficients of the helicity cross sections in eq.(|5|). The energy and momentum 
variables are defined in the text. 



Weizsacker- Williams |L5] and Altarelli-Parisi splitting functions ||16|| , 

Pih^h-x] = ^Ij^log^r (6) 

2ir 3 1 — x ra- 
ti 

As usual, x and z are the fractions of energy transferred from the t\ beam to the squark t\ 
and the gluon g after fragmentation, respectively; Q is the evolution scale of the elementary 
process, normalized by the squark mass rather than the QCD A parameter [in contrast to the 
light quark/gluon sector]. As a consequence of angular- momentum conservation, the gluon 
cannot pick up the total momentum of the squark beam. [Similar zeros have been found for 
helicity-flip fragmentation functions in QED/ QCD |Tj| |T7|.] 

By using the crossing rules {z — > 1, 1 — > x} and {1 — a; <-» 1 — x}, familiar from the 
analogous splitting functions in QED |TE[ , we derive for the elementary gluon splitting process 
into a squark-antisquark pair 

the distribution 



9 -> hh 



a 1 O 2 
P[g^t l ;x] = ^-x(l-x)\og^- (7) 



h 

after adjusting color and spin coefficients properly. This splitting function is symmetric 
under the ti <-> t\ exchange, i.e. {x <-> 1 — x}. The probability is maximal for the splitting 
into equal fractions x — 1/2 of the momenta, in contrast to spinor QED/QCD where the 
splitting into a quark- ant iquark pair is proportional to x 2 + (1 — x) 2 and hence asymmetric 
configurations are preferred. 

The above splitting functions provide the kernels for the shower expansions in perturba- 



tive QCD Monte Carlos for e + e annihilation such as Pythia [19| and Herwig 120] . They serve 
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the same purpose in the hadron-hadron versions of these generators as well as Isajet plf . 
Of course, the interpretation of the radiation processes as universal fragmentation processes 
becomes increasingly adequate with rising energy of the fragmenting squarks/gluons. 

If the ti squark is lighter than the top quark, the lifetime will be long, r > 10 _20 sec, 
since the dominant decay channel ti — > t + Xi is shut off [x? = LSP}. The decay widths 
corresponding to the 2-body decay ti — > c + Xi an d 3-body slepton decays involve the 
electro weak coupling twice and hence will be very small 0. As a result, the lifetime is much 
longer than the typical non-perturbative fragmentation time of order 1 fm [i.e. O (10 _23 sec)] 
so that the squark has got enough time to form {t\q) and {t\qq) fermionic and bosonic 
hadrons. However, the energy transfer due to the non-perturbative fragmentation, evolving 
after the early perturbative fragmentation, is very small as a result of Galilei's law of inertia. 
Describing this last step in the hadronization process of a t\ jet by the non-perturbative 
fragmentation function a la Peterson et al. |E| (which accounts very well for the heavy- 
quark analogue), we find 
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with the parameter e ~ 0.5 GeV 2 /m- . Here, x = E[{t\ q)]/E[ti] is the energy fraction trans- 
ferred from the ti parton to the (ti q) hadron etc. The resulting average non-perturbative 
energy loss 

log (- 



< 1 — x > ~ 

7T 

is numerically at the level of a few percent. 



Monte Carlo programs for the hadronization of t\ squarks link the early perturbative 
fragmentation with the subsequent non-perturbative hadronization. The relative weight of 
perturbative and non-perturbative fragmentation can be characterized by the average energy 
loss in the two consecutive steps. The overall retained average energy of the t x squarks 
factorizes into the two components, 

< x >=< x > NP < x > PT (9) 

Summing up the energy loss due to multiple gluon radiation at high energies, we find in 
analogy to heavy-quark fragmentation p3| 
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with b = (11— 2n//3)+(— 2— n//3) being the LO QCD (3 function including the colored super- 
symmetric particle spectrum. At high energies, the perturbative multi-gluon radiation has a 
bigger impact than the final non-perturbative hadronization mechanism, e.g. < x > pt k, 0.93 
for a ti beam energy E = 1 TeV and uif = 200 GeV as compared to < x > np k, 0.98. At 
low energies the two fragmentation effects are of comparable size. 



After finalizing the manuscript, we received a copy of Ref.JT0| in which the total cross sec- 



tions for squark pair production in e e~ annihilation have been discussed including squark- 
gluon and quark-gluino loops, yet not the gluon-jet distributions analysed in the present 
note. 



6 



We thank our colleagues at the LEP2 Workshop who demanded the analysis presented 
here to refine the experimental stop search techniques. 
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Figure 1: Generic diagrams for titi production in e + e collisions, (a) Born level; (b) virtual 
QCD corrections; (c) gluon radiation. 
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Figure 2: Coefficient of the QCD correction to the total cross section; shown is Pf(fi), cf. 
eq.(|), with /3=(1- 4m|/s) 1 / 2 . 
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